Polymer Effects on Heat Transport in Laminar Boundary Layer Flow 
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We consider a laminar Prandtl-Blasius boundary layer flow above a slightly heated horizontal 
plate and study the effect of polymer additives on the heat transport. We show that the action of 
the polymers can be understood as a space-dependent effective viscosity that first increases from 
the zero-shear value at the plate to a maximum then decreases exponentially back to the zero-shear 
value far away from the plate. We find that with such an effective viscosity, both the horizontal and 
vertical velocities near the plate are decreased thus leading to an increase in the friction drag and 
a decrease in the heat transport in the flow. 

PACS numbers: 



It has been known for more than 60 years that adding 
polymers into turbulent wall-bounded flows can reduce 
the friction drag significantly (see, for example, [l|, 
and references therein). A reduction in friction drag is 
equivalent to an enhancement in mass transport. Thus 
the effect of polymer additives on mass transport has 
been studied extensively during the past 60 or so years. 
On the other hand, the effect of polymers additives on 
heat transport is much less studied. Recently, an experi- 
mental study reported @ that adding polymers to turbu- 
lent Rayleigh-Benard (RB) convection of water, confined 
within a cylindrical cell heated from below and cooled on 
top, reduces the heat transport. In turbulent RB con- 
vection, there is an exact balance Q between the heat 
transport and the energy and thermal dissipation rates. 
Contributions to the energy and thermal dissipation rates 
come from the bulk of the flow as well as the boundary 
layers Q . For moderate thermal forcing as in the above- 
mentioned experimental study [using water at a moder- 
ate Rayleigh number (Ra) of 10 10 ], contributions from 
the boundary layers are dominant or significant |5fl. Fur- 
thermore, there has been experimental evidence [fj that 
the average velocity and temperature boundary layer pro- 
files in turbulent RB convection at moderate Ra could be 
described by the profiles in steady-state Prandtl-Blasius 
boundary layer flow @, H| • These experimental observa- 
tions thus suggest the possibility of understanding the 
polymer effects on heat transport in turbulent RB con- 
vection at moderate Ra by studying the effect of polymer 
additives on heat transport in laminar Prandtl-Blasius 
boundary layer flow. 

We have carried out such a study. In this Letter, we 
report our work and discuss our results. Physically, we 
can think of this boundary layer flow as the flow near 
the bottom plate of the convection cell in turbulent RB 
convection. 

For the Prandtl-Blasius boundary layer flow, the equa- 
tion of motion for the velocity field v x x + v y y is: 



Here x denotes the direction along the plate, y denotes 
the direction away from the plate, and v is the kinematic 
viscosity of the fluid. The crucial point about Eq. ([TJ is 
that the viscous term is balanced by the nonlinear advec- 
tion term and that d v 3> d x . The latter is satisfied for 



flows with large Reynolds number (Re) 
stream function 



Introducing the 



®(x,y) = VvxU (/>(£) 



with 



we obtain the famous Blasius equation 



(2) 



(3) 



Here cf>^ denotes d^cf>. The boundary conditions are v x = 
v y = at the plate and v x — > U far away from the plate, 
leading to: 

0(0) = &(0) = ; ^(oo) = 1 
Writing the temperature field as 

T(x,y)=T + (T 1 -T )8(t) 



(4) 



(5) 



where To and T\ are respectively the temperatures at the 
plate and far away from the plate (or at the center of the 
convection cell in the case of turbulent RB convection). 
Then 9 satisfies the equation 



2% + Pr 00 c = 
with the boundary conditions: 



0(0) 



1 



0(oo) = 



(6) 



(7) 



v x d x v x + Vyd y v x = vdy v x 



(1) 



and Pr = v/k is the Prandtl number. 

Here, we want to investigate the effect of polymers 
on the heat transport in this Prandtl-Blasius flow. The 
polymers produces an additional stress in the momen- 
tum equation of the fluid. This polymer stress depends 
on the amount of stretching of the polymers and is thus a 



2 



function of the dimensionless conformation tensor R\j of 
the polymers. Let the vector d denote the polymer end- 
to-end distance and pg be the polymer radius in the un- 
stretched regime, then Rij is the average over N (where 
N 1) polymers in a small region around the point 
(x, y) of the product didj/p^, i.e., = N~ 1 T,djdj / p\. 
In the simplest Oldroyd-B model of polymers J9j, the 
polymer stress is given by (v p /T)(Rij — <5y), where r is 
the (longest) relaxation time of the polymers, and v p is 
the polymer contribution to the viscosity of the solution 
at zero shear, which depends on the concentration of the 
polymers. Thus in the presence of the polymers, the 
equation of motion for the velocity field is modified by 
an additional stress that depends on R^ . By employing 
the same ideas leading to the Prandtl-Blasius equation 
([1]), we have 



v x d x v x + VydyV x = vd. 



J yy 



(8) 



Applying the transformation by £ to Eq. (JSJ does not 
generally lead to a similarity solution in that explicit ap- 
pearance of x remains in the equation for <p. This is 
known in the literature. Similarity solution have been 
obtained in some special cases with the streamwise veloc- 
ity and temperature going to some specific a;-dependent 
functions when far away from the plate [HI [H| • These 
velocity and temperature boundary conditions do not, 
however, have direct physical relevance. Here to circum- 
vent this difficulty, we recall that the Prandtl-Blasius 
flow is meant to be applicable when x is large (such that 
d y 3> d x ). Thus we make the following approximations 
by putting x = L where L 3> 1 is the length of the plate: 
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(10) 
, the 



and we have also replaced v in £ by i/q = v -+ 
zero-shear viscosity of the polymer-laden fluid. That is, 
in the presence of polymers, we have 



v x 



with polymers 



(11) 



With these approximations, the transformation of £ leads 
to a similarity solution. The resulting modified Prandtl- 
Blasius equation is: 



(1 - + 



Wi 



1 Id 



(12) 



where the Weissenberg number (Wi) and the Reynolds 
number (Re) are defined as 



tII 

Wi = — , Re = 



UL 



(13) 



and 7 = Vp/vo is a function of the polymer concentration. 
As usual in the Prandtl-Blasius approximation, all terms 
of the order of 1 /Re are nelgected in Eq. ([H 



We are interested in studying whether and how the 
heat transport is affected by the polymers. In turbulent 
RB convection, it is common to measure the heat flux 
Q in terms of the dimensionless Nusselt number (Nu), 
which is the ratio of Q to the heat flux when there is 
only conduction and is defined by 



Nu 



Q 



dT 



y=o 



2k{T x - T )/H 2(T X - T )/H 



(14) 



Here k is the thermal conductivity of the fluid, H is the 
height of the convection cell, and (■ ■ ) a is the average 
over the cross section of the cell. For the Prandtl-Blasius 
flow, taking H = L and dropping the numerical factor, 
Nu can be estimated as 



V vo 



(15) 



We want to state our major result: we will show that 
the action of the polymers is to give rise to a space- 
dependent effective viscosity that increases from the zero- 
shear value at the boundary. This increase of viscosity 
leads to an enhancement in drag, which in turn induces 
a reduction in heat flux. To solve Eq. (|T2^) . one needs to 
supplement it with specific information of R xy . In a fluid 
flow of velocity v, the dimensionless polymer end-to-end 
distance / = d/p$ obeys the differential equations 



(h — lot) + IjdjVi + thermal noise (16) 



dk _ 1 
~dt ~ 2t' 

where Iq x = cos a, Iq v = sin a and a is a random angle 
uniformly distributed in [0, 2-71"]. For the two-dimensional 
Prandtl-Blasius flow and neglecting thermal noise, we 
can rewrite Eq. (jTHJ) as 



L dl x 
U~dT 
L dl y 
U~dt 



(17) 



(Jy - lay) 
2Wi 



In order to obtain li as a function of £, we assume that 
each polymer follows the streamline of the flow such that 
d£/dt ps -{U/2L)4>. 

The quantity R xy is the average Of li over all the 
polymers contained in a small volume centered near the 
point (x,y) and over the angle a. Equations (fl"2j) . (fTTj) 
and (fTS)) are to be solved consistently. The procedure 
of how these equations are solved is rather involved and 
the details will be discussed in a forthcoming paper [l2| . 
Here we show that one can already obtain useful physical 
insights by assuming a balance, for any £, between the 
stretching of the polymers due to the velocity field and 
the relaxation, represented by the term with the factor 
1/Wi. This is equivalent to say that, in the small- Wi 
limit, we can neglect the LHS of Eqs. (JT7]) and (d3). Ne- 
glecting terms Oil/ VRe) and averaging over a, we ob- 
tain, to the order in Wi 2 : 

R xy = (1 + £<%Wi)^Wi\/Re" (19) 
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Substituting Eg. ([19]) in Eq. ([12]). we have 

2^«c + #«+2^[/(O^«] = 



(20) 



where /(£) = A£(/>££ and A = Wij is a function of Wi and 
the polymer concentration. Equation (|20[) demonstrates 
that the effect of the polymers is equivalent to introduc- 
ing a space dependent viscosity. To see this, go back to 
Eq. ([5]) and assume that v p R xy /T is equivalent to the 
product of a y-dcpendcnt effective viscosity, v v + vof(y)i 
and the velocity gradient, d y v x , we obtain: 



■ VydyV X 



vd yy v x 



-dy{[u p + vof(y)}dyV x } (21) 



which is exactly equivalent to Eq. (|20) once the rescaling 
in the variable £ is performed. 

We solve Eq. for A = 1, and show /(£) in Fig.Q] It 
can be seen that /(£) increases from zero at the boundary 
and then decreases almost exponentially towards zero. 
This increase of the effective viscosity near the plate al- 
lows us to predict the effect of the polymer on the flow. 
First, we expect an increase in the effective viscosity near 
the plate would give rise to an increase of the drag. We 
measure the drag by the drag coefficient: 



C 



v d y v x 



1 



U 2 /2 



2VRc 



(22) 



and compare C to the value Co for the flow without poly- 
mers. We expect C/Co would be greater than 1 for all 
A > 0. Such an enhancement in drag would in turn 
result in a reduction of the mass throughput in the hor- 
izontal direction and thus a decrease in the horizontal 
velocity near the boundary. These effects are indeed ob- 
served. In Fig. Q] we show the horizontal velocity v x and 
respectively for the case without polymers (denoted 
by the superscript 0) and with polymers (denoted by the 
superscript p), and also compare the vertical velocities, 
vP and Vy, with and without polymers. It can be seen 
that in the presence of polymers, both the horizontal and 
vertical velocities decrease in the region near the plate. 

Next, we show that drag enhancement implies a reduc- 
tion of the heat transport or Nu. Upon double integration 
of Eq. ((6]) by £, we obtain: 



A(0) 



ds\ exp[— Pr 



fas 2 )ds 2 /2] 



(23) 



which tells us that Nu is a functional of $, where $(£) = 
$ <}>(s)ds. We can thus calculate JNu, the variation in 
Nu induced by a variation in $, <5<& (^) : 



(5Nu 



Pr Nu 2 



cxp[-Pr $(s)/2] <5$(s) ds (24) 



Since v x = Ufa, the mass throughput in the x direction 
across a distance £ is given by Ufa£), and thus a reduc- 
tion in mass throughput implies an <5<I> < 0. Therefore, 
Eq. (|24p shows that drag enhancement implies a reduc- 
tion in Nu or the heat transport. 




FIG. 1: Viscosity and velocity profiles obtained by solving 
Eq. p0[) with A = and A = 1. We show the space-dependent 
viscosity /(£) for A = 1 (squares), the horizontal velocity, 
v x(0/ u (triangles) and v%(£)/U (circles), without (A = 0) 
and with polymers (A = 1), and the difference in vertical 
velocity, [v*(£)-v°(Z)]/(U/2\/Re) (solid line). 



The above considerations are not restricted to the par- 
ticular form of /(£) = A^cf>^ obtained in the small- Wi 
limit but are true for a generic form of /(£) which displays 
the same features of first increasing near the boundary 
to some maximum then decreasing to zero far away from 
the boundary as shown in Fig. [TJ In particular, we show 
that this is the case by using a different form of /(£): 



/(£)=A£exp(-0 



(25) 



that has similar qualitative features. The constant A 
relates the polymer characteristics to the conformation 
tensor Rij , and generally increases with Wi and the con- 
centration of the polymers. 

Solving Eq. (|2"5j) with Eq. (|2"fj|) . we again find an in- 
crease in drag and a reduction in heat transport. More- 
over as shown in Fig. [2l the extent of the effect is in- 
creased by increasing A, i.e., by increasing the effect of 
the polymers into the system. For the range of A shown 
in Fig. [3J the amount of the reduction of heat flux is 
quantitatively consistent with that observed in Ref. Q. 
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FIG. 2: Dependence of C '/Co, (triangles) and Nu/Nuo, (cir- 
cles) on A for /(£) given by Eq. 1(25)1. 
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FIG. 3: jg (£) as a function of Wi at fixed 7 = 0.2 obtained by 
the numerical solution of the full nonlinear problem, Eqs. (|12|) , 
(|17|l and (|18|l . Wi increases from bottom to top, taking the 
values 0.5, 1.0, 1.5, 2.0, 2.5, and 2.8 respectively. 

In the remaining of this Letter, we show that the 
results that we have discussed capture the qualitative 
physical picture of the full nonlinear problem, defined 
by Eqs. dHJ) , dTTJ) and (H]). We have solved the full 
problem, by integratin g nu merically the equations using 
an iterated procedure [12j. In particular, we show that 
Rxy = [1 + WiVRe for some function g{£). Sub- 

stituting this into Eq. (fT2"]) gives 

2&« + #« + 2d ( [ 7 ff = (26) 

Comparing Eq. ((26)) with Eq. ([20]) for the small- Wi limit, 
we see that 7ff(£)] plays the role of /(£). In Fig. 02 we 
show 7<?(£) for different values Wi at a fixed value of 7 = 
0.2. Notice that the functional shape of <?(£) is similar 
to the one used in Eq. ([25]) . Not surprisingly, we find 
drag enhancement and heat transport reduction for the 
full nonlinear problem as shown in Fig. @] 




Wi 

FIG. 4: Drag enhancement and heat transport reduc- 
tion for the full nonlinear problem and the dependence of 
C/Co (squares) and Nu/Nuo (circles) on Wi. 

In summary, we have shown the following results: 

(1) It is possible to write a self consistent set of equa- 
tions for Prandtl-Blasius boundary layer flow with 
polymers. 

(2) The physical effect of the polymers is equivalent 
to introducing a space-dependent effective viscosity 
that is increasing near the boundary. 

(3) An increasing viscosity near the boundary leads to 
drag enhancement, which in turn induces a reduc- 
tion in heat transport. 

Our theory may thus explain the recent experimental 
results observed in turbulent RB convection [3| . 
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